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Water wave propagation and scattering over topographical bottoms 
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Here I present a general formulation of water wave propagation and scattering over topographical 
bottoms. A simple equation is found and is compared with existing theories. As an application, the 
theory is extended to the case of water waves in a column with many cylindrical steps. 



PACS numbers: 47.10.+g, 47.11. +j, 47.35. +i 



I. INTRODUCTION 

There have been many approaches for investigating 
propagation of water waves over various bottom topogra- 
phies. A vast body of hterature exists. For brevity, I refer 
the reader to the textbook ^ . Here I would hke to derive 
from the first principle a simple but coherent formulation 
for the problem. It will be shown that this approximate 
approach compares favorably with existing approxima- 
tions when applied to the cases considered previously. 
The advantage of the present approach is obvious: it is 
simple, accommodating, systematic, and can be easily 
numerically programmed. In particular, here I explicitly 
show that it respectively recovers three previous results 
for shallow water, deep water, and scattering by rigid 
cylinders standing in water. I will first give a theory 
for general bottom topographies. Then I will extend to 
study the case of water wave propagation and scattering 
in a column with many cylindrical steps. 



II. GENERAL THEORY 

Consider a water column with an arbitrary bottom to- 
pography. We set up the coordinates as follows. Let the 
z axis be vertical and directed upward. The x — y plane 
rests at the water surface when it is calm. The depth of 
the bottom, which describes the bottom topography, is 
denoted by h{x,y), and the vertical displacement of the 
water surface is rj(x,y,t). Now we derive the governing 
equations for the water waves. 

Consider a vertical column with a base differential el- 
ement dxdy at (x, y). The change rate of the the volume 
of the column is 



d_ 

dt 



ri{x, y, t)dxdy. 



By conservation of mass, this would equal to the net vol- 
ume flux from all the horizontal directions, i. e. 



d_ 

dt 



r]{x,y,t)dxdy 



dzvj_{x, y, z, t) 



where = {dx,dy), and '_L' denotes the horizontal di- 
rections. This gives us the first equation 



9 , 

-g^mx,y,t) 



dzvj_{x, y, z, t) 



(1) 



The second equation is obtained from the Newton's 
second law. From the Euler equation for incompressible 
ideal flows 

dtv + [v ■ \7)v = — Vp — pgz, 
P 

which is valid at z = 0, with g being the gravity acceler- 
ation, and 



we obtain 
d 



—v^_{x,y,0,t) + [{v-V)v] 



(2) 



=0 = -9^±V- (3) 



Note when the liquid surface tension is included, the fol- 
lowing term should be added to Eq. (0) 



(4) 



in which cr is the surface tension coefficient. In this paper, 
for short, we ignore this effect. 

Another equation is from the boundary condition at 
z = h, which states 



\z— — hi^x,y) 



0, 



(5) 



where n is a normal to the bottom. For an incompressible 
fluid, wc also have the following Laplace equation. 



V • v{x,y,z,t) = 0, 



(6) 



in the water column. 

Equations (|l|) , (||) , (^) , and (j^) are the four fundamen- 
tal equations for water waves. 



Linearization 



For small amplitude waves, i. e. j] << h, we can ignore 
dxdy, f^ijQ non- linear terms in (|^) and (^. Such a linearization 
leads to the following two equations 
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—r]{x,y,t) = -Vj. 



h{x,y) 



dzv±{x, y, z, t) 



(7) 
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and 



-g^v±{x,y,0,t) + gV±r]{x, y, t) 



0. 



(8) 



These two equations together with Eqs. (g) and de- 
termines scattering of water waves with a bottom topog- 
raphy. 



B. Propagation approximation 

Here we provide an approximate solution to Eqs. (|^), 
(|), (|), and (|). The procedure is as foUows. When 
the variation of the bottom topography is smaller than 
the wavelength (to be determined self-consistently) , we 
can first ignore terms involving Vj^/i, and solve for the 
velocity field. For the incompressible fluid, the velocity 
field can be represented by a scalar field, i. e. 

v{x,y,z,t) = W^{x,y,z,t). 

We write all dynamical variables with a time dependence 




Transmission 
Reflection 



FIG. 1; Transmission and reflection coefficients versus (3 — kh 
for an infinite step, obtained from Eq. (^). While the result 
for the reflection agrees very well with that in Refs. 
there is some discrepancy in the transmission results within 
the range of kh between 0.4 to 1.2; the largest discrepancy of 
about 15% occurs around kh = 0.8 for the transmission. The 
legends are adopted from H 



(this time fact is dropped afterwards for conve- Here </> and are determined by 



nience). This procedure leads to the following equations 
for $. 



(vi + P)0 = o, 



with 



VHix,y,z)^0. 



d 



cc;^$(f,0)-f.9 — $(f,0) = 0; (f=(a;,y)). 



(9) 



(10) 



and 



(Vi 



(16) 



(17) 



The first approximation is made at the bottom (z — —h). 
The boundary condition at the bottom reads 

— $(f,-/i) = — $(f,-/i) + Vi-Vi$(f,-/i) = 0. (11) 

We approximate that fi is in the z direction by neglecting 
the second term in the above equation. This is valid as 
long as Vj_/i << kh. Thus the bomidary condition gives 



d 



(12) 



Note that this condition is exact in the case of step-wise 
topographical bottoms, to be discussed later. Eqs. (||), 
(|l^), and (|l|) lead to the solution for $ 

^{x,y,z) = (l){x,y) cosh{k{z + h)) + 

^<j)n{x,y)cos{kniz + h)), (13) 

n 

where k satisfies 

= gk{x, y) tanh(fc(x, y)h{x, y)), (14) 
and kn satisfies 



Eq. ( p7| ) leads to evanescent wave solutions. 

The second approximation is to ignore the summation 
terms in Eq. (p^. Such an approximation is based upon 
the following consideration. The summation terms repre- 
sent the correction of evanescent waves caused by irregu- 
larities such as sudden changes of depth. As these waves 
are spatially confined, it is reasonable to expect that such 
a correction will not affect the overall wave propagation, 
and the general features of the wave propagation. Indeed, 
when we apply the later approximate solution to the ex- 
treme case of propagation of water waves over an infinite 
step, we find that our results agree reasonably well with 
that from two other approximate approaches^, |j. For 
example, the difference in the reflection results is uni- 
formly less than a few percent for a wide range of fre- 
quencies. The largest discrepancy can happen for the 
transmission results, but the difference is still less than 
15%. Furthermore, we find that the derived result is 
in agreement with that of Kirby for the case of waves 
over a flat bed with small ripples As matter of fact, 
in this case, it can be shown that after a mathematical 
manipulation 0] Eq. (2.11) in becomes essentially the 
same as the following Eq. (pl|). 

Under the above approximations, we have 



^{x, y, z) ~ (j)(^x, y) cosh(A;(z -I- /i)), 



(18) 



and 



= gkn{x, y) tan(fc„(a;, y)h{x, y)). 



(15) 



Vi_(x, y, z) w cosh(fc(z + /i))Vj_0. 



(19) 
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Now taking Eqs. ^ and (|l9|) into Eqs. (|) and (|), we 

get 



D. The situation of deep water or high frequencies 



(20) 



For convenience, hereafter we write as V when it acts 
on the surface wave field rj. That is 



V ( ^Vry ) + = 0, 



where k satisfies 

uj'^ = gk{f)ta,nh{k{r)h{f)). 



(21) 



(22) 



(23) 



From this equation, we can have the conditions hnking 
domains with different depths as follows: both rj and 
tanh(fc/i) ^ _ i^jy g^j-g continuous across the boundary. 

Eq. (|2^) is similar to what is known as the mild-slop 
approximation |l| : 



iv(i,v,)+, = o, 

where c is given by 

1 



2kh \ 
smh{2kh) ) 



(24) 



(25) 



Eq. (|2J) was derived by a number of authors under the 
situation that V/i << kh. In fact, under this condition it 
can be shown that Eq. (^ and Eq. (|2^) are equivalent. 

Note that when the surface tension is added, Eq. ( |2l| ) 
becomes 



tanh(fcft,) 



V 77 



gp 



H r/ = 0, (26) 



with Eq. (p3[) becoming 



cj^ = [gk + -k^ ) tanh(fc/i). (27) 



C. The situation of shallow water or low 
frequencies 

In the case of shallow water, i. e. fc/i << 1, we obtain 



from Eq. (21) 



For the deep water case, hk » 1, we have 

k=—, 

9 

and 

,4 
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(29) 



(30) 



In the deep water, the dispersion relation is not affected 
by the bottom topography. 



E. Scattering by infinite rigid cylinders 

Equations (|2^) or ( p^ are also applicable to another 
class of situation which has been widely studied in the 
literature. That is, the scattering of water waves by infi- 
nite rigid cylinders situated in a uniform water column. 
When applying ( pl| ) or ( p2| ) to this case, we find that 
these two equations are actually exact. In the medium, 
the wave equation is 



(V + k^)T] ^ 



(31) 



with the boundary condition at the i-th cylinder 

■ V77I, = 0, (32) 

obtained as we set the depths of the cylinders equal zero; 
fii is a normal to the interface. In fact, in this case, 
the problem becomes equivalent to that of acoustic scat- 
tering by rigid cylinders, and all the previous acoustic 
results will follow ^, |l^, such as the interesting phe- 
nomenon of deaf bands. 



III. WATER WAVES IN A WATER COLUMN 
WITH CYLINDRICAL STEPS 

The problem we are now going to consider is illustrated 
by Fig. 1^. We consider a water column with a uniform 
depth h. There are N cylindrical steps (or holes when 
hi > h) located in the water. The depths of the steps 
are measured from the water surface and are denoted 
by hi and the radii are a^. In the realm of the linear 
wave theory, we study the water wave propagation and 
scattering by these steps. 



Band structure calculation 



V • (Wr;) + — 7] = 0. 
9 



(28) 



This is the fundamental equation governing the small am- 
plitude waves in shallow water, first derived by Lamb[^. 



When all the steps are with same hi — h2 — ■ ■ ■ — h^ 
and the radius a, and are located periodically on the bot- 
tom, then we can use Bloch's theorem to study the water 
wave propagation. Assume the steps are arranged either 



4 



and is the structure factor 




FIG. 2: Conceptual layout (side view of the three dimensional 
coordinates) : There are A'' cylindrical steps located in a water 
column with depth h. The depths of the steps are denoted 
hy hi{i = 1,2, ... , N) measured from the upper surface of the 
water column, and the radii of the steps are denoted by a^. 
The coordinates are set up as shown. The steps are located 
at fi. The y-axis lies perpendicularly to the page 



in the square or hexagonal lattices, with lattice constant 
d. Here we use the standard plane- wave approach |llj p^ . 
By Bloch's theorem, we can express the field 77 in the fol- 
lowing form 



(33) 



where r = {x, y), G is the vector in the reciprocal lattice, 
and K the Bloch vector. 

In the present setup, the bottom topograph is periodic, 
so we have the following expansion 



tanh(fcft,) 
k 



(34) 



with 



/tanh(fci/ii) tanh(fcft,)\ tanh(fc/i) 
k, + 



for 



and 



G = 0; 



^ ^ ^ / tanh(fcifei) _ tanh(fcM\ ^ 
\ ki k J 



for 



Here ki and k are determined by 

up' = gfci tanh(fcift,i) = gfctanh(/c/i), (37) 
and fs is the filling factor given by[p^ 



Square lattice 



I ■^ "(^)^' Hexagonal lattice, 



F,{G)^2f, 



\G\a 



Substituting Eqs. (||) and (|3|) into Eq. (|T|), we get 
^gg<g,(i?,c^)C(G',i?) = 0, (38) 



G' 



with 



Qq q. {K, c.) = [(G + K) ■ (G' + K)\A{G - G') -—^qq, 
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The dispersion relation connecting K and u is deter- 
mined by the secular equation 



dot 



[G + k)- (G' + K)]A{G - G') - —h,G' 



0. 



G,G' 



(39) 

For the shallow water, we have tanh(fcft,) sa kh, and 
thus tanh(fc/i)/fc w h, then by 



with 



AiG) 



h{x,y) = Y^A{G)e^^-\ 



{hi - h)f, + h, for G = 0; 
{hi^h)F,{G), for G 7^0. 



(40) 



(41) 



B. Multiple scattering theory 

The shallow water wave propagation in the water col- 
umn with cylindrical steps can also be investigated by 
the multiple scattering theory. Without requiring that 
all the steps are the same, we can develop a general for- 
mulism. 

In the water column, the wave equation reads 

(V^ + fc2),7 = 0, (42) 
with k being given by 

= gktanh{kh) 
Within the range of the i-th step, the wave equation is 
(V^ + fcf),y, = 0, (43) 

with 

uj"^ — gkta.nh{kh) 
At the boundary of the step, the boundary conditions are 



ta,nh(kihi) 



n ■ Vrii 



tanh(A;i hi) 



h ■ V77 



(44) 
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derived from the conservation of mass, and 



(45) 



Here dfli denotes the boundary, and n is the outward 
normal at the boundary. 

Equations M2|) and (|4^) with the boundary conditions 
in (44) and ( [45D completely determine the shallow wa- 
ter wave scattering by an ensemble of cylindrical steps 
located vertically in the uniform water column. By in- 
specting, we see that this set of equations is essentially 
the same as the two dimensional acoustic scattering by 
an array of parallel cylinders We following |g| to 

study the scattering of shallow water waves in the present 
system. 

Consider a line source located at fg. Without the cylin- 
der steps, the wave is governed by 



(V^ + fc2)G(r - rs) = -47r5(2)(f - f,), 



(46) 



where H^^^ is the zero-th order Hankel function of the 
first kind. In the cylindrical coordinates, the solution is 



(47) 



In this section, 'i' stands for 

With N cylinder steps located at n {i — 1, 2, • • • , iV), 
the scattered wave from the j-th step can be written as 



^ i7rAiHi'\k\r-r,\)e'-f^^ 

n— — oo 



(48) 



where Hn ^ is the n-th order Hankel function of the first 
kind. is the coefficient to be determined, and (j)r-r- 
is the azimuthal angle of the vector r — fi relative to the 
positive X-axis. 

The total wave incident around the i-th scattcrcr 77*^ (r) 
is a superposition of the direct contribution from the 
source rjo^f) = G{f — rs) and the scattered waves from 
all other scatterers: 



N 



Vlnir) ^ Voir) + ^^^^ 



(49) 



In order to separate the governing equations into modes, 
we can express the total incident wave in term of the 
modes about rv: 



(50) 



The expansion is in terms of Bessel functions of the first 
kind J„ to ensure that ?7i„(r) does not diverge as r ^ r^. 
The coefficients are related to the in equation 
( ^ ) through equation (^9[). A particular represents 
the strength of the n-th mode of the total incident wave 
on the z-th scatterer with respect to the i-th scatterer's 
coordinate system (i.e. around fi). In order to isolate 



this mode on the right hand side of equation (^), and 
thus determine a particular i?^ in terms of the set of A^, 
we need to express 77s (r,rj), for each j 7^ i, in terms of 
the modes with respect to the i-th scatterer. In other 
words, we want rjsir^fj) in the form 



,(f,r-)= J2 C^/Jn(.k\r-n\)e 



l<Pr 



(51) 



This can be achieved (i.e. C^'' expressed in terms of 
through the following addition theorem[|4|: 



Ui{k\r-r-l\)e''^^^-'-^.{52) 



1 — — OC 



Taking equation ( p2| ) into equation (p8[), we have 

00 



E HiZik\rl-r,\)e 



^ J,(fc|f-rl|)e"*^^-'-^..(53) 



Comparing with equation we see that 

00 

Ci^^= E i^A^^l^(fcK-r}|)ei('-")*--^ (54) 



Now we can relate Bl^ to C^'^ (and thus to Aj) through 
equation (^9|). First note that through the addition the- 
orem the source wave can be written. 



77o(r) = i7riJ^^^(fc|f- r; 



where 



Er=-^^SlMk\r-n\)e 



Si^mH[^Xk\rl-rs\)e-''^^--^ 



(55) 



(56) 



Matching coefficients in equation (49) and using equa- 
tions (|50|), ( |5l| ) and (55), we have 



N 



— E 



(57) 



or, expanding C,^'*, 



N 00 

(58) 

At this stage, both the 5*,^ are known, but both i?^ and 
Aj are unknown. Boundary conditions will give another 
equation relating them. 
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The wave inside the i-th scatterer can be expressed as 

oo 

Vlnti^^ DUnihlr-nDM^-^-l. (59) 



n—~oo 



Taking Eq s. (|4^ ), js^) and ( |59| ) into the boundary con- 
ditions in ([44|) and ([t5[), we have 



(60) 



(61) 

where "' refers to the derivative. Ehmination of D^^ gives 



(62) 



where 



tanh(/c/i) TjW'/ 



- Ts e 



-in0j,-. 



tanii(fci;ii) J'n{kai)Jn{kiai) Jn{kai)J!^{kiai) 
If we define 

= 5;/i^ = ijiV. 

and 

then equation ( |5^ ) becomes 

Af oo 



(63) 
(64) 

(65) 
(66) 



If the value of n is hmitcd to sonic finite range, then this 
is a matrix equation for the coefhcients A]^. Once solved, 



the total wave at any point outside all cylinder steps is 



?7(f) = i^i/W(fc|f-r=;|) + 

N oo 

^ J2 i^A;ijW(fc|r-rl|)ei"'^--. (67) 

i—l n— — oo 



We must stress that total wave expressed by eq. (|67| ) 
incorporate all orders of multiple scattering. We also 
emphasize that the above derivation is valid for any con- 
figuration of the cylinder steps. In other words, eq. ( |6^ ) 
works for situations that the steps can be placed either 
randomly or orderly. 

For the special case of shallow water {kh << 1), we 
need just replace in Eq. (|6^) by 



Hl^\ka,)Jl^{ha,) - ^-^Hk'>'ika,)J„{ha,) 
h^Jn(kai)Jn{kiai) - Jn{kai)J^{kiai) 



F' = 



(68) 



IV. SUMMARY 



In summary, here we have presented a general the- 
ory for studying gravity waves over bottom topographies. 
The results have been extended to the case of step-wise 
bottom structures. The model presented here is simple 
and may facilitate the research on many unusual wave 
phenomena such as wave localization flS, n&. 
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